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<^ Abstract 

Analytical formulae for the points and weights of two fifth-order 
i -^H quadrature rules for C3, the 3-cube, are given. The rules, originally 

formulated by A. H. Stroud in 1967, are discussed in greater detail 
^ in terms of both the setup of the basic equations and the method of 

obtaining their solutions analytically. The primary purpose of this 
T— I paper is to better document what we feel is a particularly practical 

^ quadrature rule (e.g. in finite element calculations) and one for which 

we felt comprehensive information was scarce. 
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§ 1 Introduction 

• ^ In 1967, A. H. Stroud published an article [8] on fifth-degree integration 

^ formulas for several symmetric, n-dimensional regions. In the first sentence 

of Section 2 of that article, Stroud mentions that "Unless stated otherwise 
we assume that n > 4." He goes on to give a general description of what 
is now a well-known method for determining non-product quadrature rules 
for several standard regions including the n-cube, n-sphere, and the entire 
n-space. 

This description is followed by a number of tabulated quadrature rules 
for specific n having about six digits of precision. A particularly interesting 
rule (herein referred to as "Stroud's first rule") is given for the 3-cube, "C3". 
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In the last sentence of the paper, Stroud states: "Previously no such 13 point 
formula was known for C3." A second rule (herein referred to as "Stroud's 
second rule") for the 3-cube, having some points outside the region, is given 
as well. 

Unfortunately, the generic equations for the ra-cube given by Stroud in [8] 
do not apply to the n = 3 case. This fact makes it impossible to reproduce 
Stroud's tabulated results (in order to e.g. compute the points and weights 
to a higher precision) without investing some amount of time redoing the 
algebra oneself for the n = 3 case. In Stroud's famous 1971 compendium on 
quadrature rules [9] , the first rule for the 3-cube is given in additional detail: 
fifteen decimal digits of precision (suitable for double-precision calculations) 
are given for the points and weights, and an eighth-order polynomial is given 
for determining two of the rule parameters, but no additional guidance is 
given in this short table entry. Stroud's second rule for the 3-cube is also not 
given in detail in [9]. 

In this paper we provide additional details on Stroud's first and sec- 
ond rules for the 3-cube. In particular, we give analytical formulae for the 
points and weights, which, to our knowledge, have not been previously pub- 
lished. Although we came across a number of potential references for these 
rules [1, 2, 3, 4, 5, 7], none were found that contained exactly the information 
given here. The purpose of this paper is therefore to provide an accessible 
reference for what we feel is an especially practical rule. For example, the 
analytical formulae presented here are now being used in the general purpose 
finite element library, LibMesh [6]. The reasons for preferring an analytical 
solution over tabulated numerical values are obvious, chief among them be- 
ing adequate precision in any computing environment, even those not yet in 
existence. 



2 Basic Equations for Stroud's Fifth-Order 
Rules for the 3-cube 



Stroud's technique for obtaining fifth-order quadrature rules for symmetric 
regions R is summarized as follows: choose points z/j and weights Ai such 
that the approximation 




(1) 



2 



1 2 * 



for which 



is exact for all monomial functions / := 

|a| := «! + ^2 + • • • + «n < 5 (2) 

Due to the symmetry of the regions studied, Stroud additionally restricts the 
set of possible rules to only those which contain symmetric pairs of points 
ii/j and their negatives, both having weight Ai. (If i/j is the origin, then its 
negative is not included in the rule.) The rules are also restricted to have 
M = ^{n^ + n) + 1 distinct points. Such rules are then automatically exact 
for all monomials in which |a| is odd. 

The condition that Eqn. (1) be exact for monomials with \a\ = 0,2,4 is 
then equivalent to the matrix equation 

X^AX = C (3) 

which is Stroud's Eqn. (5), where A := diag^Ai, A2, Am), X contains 
various (quadratic) products of the z/,, and C contains corresponding exact 
monomial integrals for the given region. (The reader should refer to [8] for 
additional details.) Assuming C is non-singular, we can take inverses and 
thus Eqn. (3) is equivalent to 

XC-^X^ = A-^ (4) 

which is Stroud's Eqn. (7). Finally, Stroud chooses the M points and weights 
in a special way by taking (for the particular case n = 3, M = 7) 

= {v,v,v) 

^2 = (A,e,0 1 

^^3 = (^, A, > weight B vq = (/i, 7, ^) ) weight C 
^4 = (e,e,A) J 

With this special choice of points, and taking := 0, Eqn. (4) has the left- 
hand side 
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(5) 



and the right-hand side 



A-^ = 32 diag {A-\ B-\ B-\ B-\ C-\ C'^) 



(6) 



Wc have added extra horizontal and vertical lines to the matrix given in 
Eqn. (5) to emphasize the symmetric sub-structure present in the govern- 
ing equations. In Eqns. (5) and (6) we have also introduced the following 
submatrices: 



Mi{x,y) :-- 



mi 
mi 
mi 



mi{x,y) := 19 - ISx^ - 30y^ 



1712 n^3 n^3 
rris m-2 rn^, 
ms ma m2 



(7) 



m2{x, y) 45x^ - SOx^ - GOy^ + 126y^ + 72x'^y'^ + 19 (8) 
m3(x, y) := 45y^ - SOx^ - GOy^ + 126x^y2 + 72xy^ + 19 (9) 



m4 m4 m5 
m4 m5 m4 
ms m4 m4 



m4 : 



m5 : 



= 45{V^ - 30C^ + 45/^2 - 15A=^ + 457^^^ + SGA^/x^ 

+ 36^V7 + 36A^^7 - SO/x^ - 157^ + 19 (10) 
126CV^ - 30^2 + 457^A2 - ISA^ + 72A^/X7 
- 30/^2 - 157^ + 19 (11) 



The special form assumed for the quadrature points and weights has 
reduced the size of the original system of equations quite drastically, from 
7^ = 49 equations to fewer than 10. 
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3 Analytical Solution of the Equations 



The equations defined by Mi and Af 2 may be used to determine all possible 
solutions for the triplets (A, ^, B) and (7, /x, C) independently and simulta- 
neously. The following procedure is used: we rearrange Eqn. (7) to solve for 
y (resp. x) and insert it into the rris equation, (9). Since the equation for rris 
has odd powers of both x and y, wc obtain two possible forms of the ma equa- 
tion: one for the positive square root of x (resp. y) and one for the negative 
square root. Setting the product of the positive- and negative- root versions 
of Eqn. (9) equal to zero leads to an eighth-order polynomial equation in x 
(resp. y) which is equivalent to a quartic polynomial equation in x"^ (resp. y'^). 
Wc can solve this quartic equation for x"^ (resp. y^) analytically, and, finally, 
the resulting (a;^,7/^) pairs may be substituted into the 1712 equation, (8) to 
solve for the weights. 

The eighth-order polynomial in x (quartic in x^) arising from the previously- 
described procedure is given by 

o 4 2828796 , , , 

1330425x^ - 3108780x^ + 2339622a;^ x^ + 361 = (12) 

5 

while the eighth-order equation for y is 

„oi7 8 363204 833454 , 30324 

53217y" — y° H ^^y —y + 361 = (13) 

We can use any suitable CAS to solve Eqn. (12) and obtain the four 
solutions of (recall that our generic x variable corresponds to either A or 
7) analytically. The result is 

2 _ r 2 2 1 

1919 + 148^ ± 4i_ 1919 - 148v^ ± 4i+ , 

^ (14) 



3285 ' 3285 

where the short-hand notation 



y 71440 ± 6802v^ (15) 

is used, and will also be used throughout this paper. In Eqn. (14) and those 
which follow, for any quantity with an i,j subscript and an ambiguous "±" 
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sign, the i subscript always refers to the top sign throughout the equation, 
while the j subscript refers to the bottom sign. 

In a similar manner, we obtain the following analytical solutions for y'^ 
(recall that our generic y variable corresponds to either ^ or /i) by solving 
Eqn. (13) 

y"^ = {z/L. 2/3,4} 

_ J112I -74v^^2t_ 1121 + 74v^^2t 
~ I 3285 ' 3285 

(Note: we have ordered the and yf solutions such that pairs (xj, ?/«) satisfy 
the nil equation (7).) We can now substitute the {xi,yi) pairs given above 
into the generic m2 equation (repeated here) 

45x^ - 30x^ - 60y^ + l2Qy^ + 72x^y'^ + 19 = 32^-^ (17) 

to solve for the generic weights w. The weights Wi (which we have already 
scaled by |, since each weight solved for in Eqn. (17) is actually twice the 
true value) correspond to the B and C parameters in the original equations, 
and are given by 

w = {^1,2,^3,4} 
_ J 133225 
~ [ 260072 + 1520v^ ± (l33 + 37v^) t_ ' 

133225 

260072 - 1520v^ ± (l33 - 37v^) t 

Upon simplification, the equation, (11) yields 

126^V^ + 4572A2 + 72A^/i7 + 19 = (19) 

This implies that exactly one of the four parameters A, ^, /z, 7 must be 
negative. (Note: Conversely, three of the four parameters could instead be 
negative, but since the rule always includes a point and its negative, this 
is equivalent to one of the four parameters being negative.) This means 
that, when taking square roots of the x^ and y"^ values obtained previously, 
exactly one negative root must be selected. In both of his rules, Stroud [8] has 
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selected ^ as the negative root, and we shall follow the same convention here. 
Using the rest of Stroud's tabulated results as a guide, we have compiled 
Table 1, which gives the corresponding analytically-obtained Xi, yi, and Wi 
values. It may of course be verified that these analytical values also satisfy 
Eqns. (10) and (11). 

Table 1: The middle column gives numerical approximations to the analytical 
results obtained in this section, for reference. In the right-hand column, we 
give the corresponding parameter originally obtained by Stroud [8]. The 1 
and 2 subscripts in the third column correspond to the first (with all points 
inside the region) and second (with some points outside the region) rules 
reported by Stroud. 



Xi 


1.0146309695 


72 


X2 


0.7291297984 


A2 


Xs 


0.8803044067 


Ai 


X4 


0.0252937117 


7i 


yi 


0.3443767286 


/^2 


y2 


0.6062327951 


-6 


y3 


0.4958481715 


-6 


Z/4 


0.7956214222 


/il 


Wi 


0.4075948702 


C2 


W2 


0.6450367090 


B2 


W3 


0.5449873514 


Bi 




0.5076442277 





4 Summary 

To summarize the analytical results reported, we give both numerical approx- 
imations and analytical forms for the various parameters in this section. In 
Tables 2 and 3 the numerical and analytical values, respectively, for Stroud's 
first quadrature rule are given. The numerical approximations are given to 32 
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digits of accuracy for convenience and because in the case of the second rule, 
such highly-accurate values have not been previously published. Stroud's 
second fifth-order rule (which is less useful for finite element calculations due 
to the fact that some of the points lie outside the region of integration) is 
likewise summarized in Tables 4 and 5. 

Table 2: Values for Stroud's first fifth order quadrature rule for the 3-cube, 
originally reported in [8], to 32 decimal digits. 

7] O.OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOE+00 

A 8.80304406699309780477378182098603E-01 

^ -4.95848171425711152814212423642879E-01 

/i 7.95621422164095415429824825675787E-01 

7 2.52937117448425813473892559293236E-02 

A 1.68421052631578947368421052631579E+00 
B 5.44987351277576716846907821808944E-01 
C 5.07644227669791704205723757138424E - 1 
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Table 3: Analytical representations for Stroud's first fifth order quadrature 
rule for the 3-cube, originally reported in [8]. The results correspond to the 
numerical approximations given in Table 2. When taking the square root, 
the positive root is always assumed unless specified otherwise. 



A 



7 



1919 


- 148 




+ 4t+ 




32 


35 




1121 


+ 74v 


/l9 - 


2t+ 




328 


5 




1121 


+ 74v 


/l94 


-2t+ 




328 


5 




1919 


- 148 


Vl9 


-4t+ 




32 


35 





A 



B 



C 



32 
19 



133225 



260072 - 1520v^ + (l33 - 37v^) U 
133225 



260072 - 1520v^ - (l33 - 37v^) t+ 



t± := V 71440 ±6802V19 
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Table 4: More accurate values for Stroud's second fifth order quadrature rule 
for the 3-cube, originally reported in [8], to 32 decimal digits. 



7] O.OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOE+00 

A 7.29129798350178619517304350180274E-01 

^ -6.06232795147414690867819787871004E-01 

fi 3.44376728634554308789703229927940E-01 

7 1.01463096947441524044348381753366E+00 

A 1.68421052631578947368421052631579E+00 

B 6.45036708927015064146303620658267E-01 

C 4.07594870020353356906327958289101E-01 
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Table 5: Analytical representations for Stroud's second fifth-order quadrature 
rule for the 3-cube, originally reported in [8]. The results correspond to the 
numerical approximations given in Table 4. When taking the square root, 
the positive root is always assumed unless specified otherwise. 



A 



7 



1919 + 148 




-4t_ 


32 


B5 




1121 - 74v 


/T9 + 


2t_ 


328 


5 




1121 - 74v 


/T9 - 


2t_ 


328 


5 




1919 + 148 


Vl9 


+ 4t_ 


32< 


35 





A 



B 



C 



32 
19 



133225 



260072 + 1520v^ - (l33 + 37v^) t_ 
133225 

260072 + 1520v^ + (l33 + 37v^) t_ 



t± := V 71440 ±6802V19 
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